УТОЧНЕННАЯ МЕТОДИКА РАСЧЕТА ОПТИМАЛЬНОГО ДИАМЕТРА ГАЗОПРОВОДА 

С. В. Голдаев 

Томский политехнический университет 
Тел.: (382-2)-5б-40-10 

Предложена уточненная методика расчета оптимального диаметра газопровода, в которой учтена зависимость коэффициента 
гидравлического сопротивления от диаметра трубы и режима движения газа. Показано, что приближенная методика, предпо- 
лагающая постоянным упомянутый коэффициент, дает на 10... 12 % завышенные значения оптимального диаметра и на 2...3 % — 
суммарные затраты на строительство и эксплуатацию газопровода. 


Доставка и распределение природного газа пот- 
ребителям (промышленным, энергетическим, бы- 
товым) требуют больших затрат денежных и мате- 
риальных средств, в том числе одного из наиболее 
дефицитных видов проката - труб. Поэтому сни- 
жение стоимости и металлоемкости систем газос- 
набжения имеет важное значение [1]. При транс- 
портировке газа по трубам преодолеваются гидрав- 
лические сопротивления, на что тратятся опреде- 
ленные энергетические затраты. Кроме того, под- 
держание газопровода в рабочем состоянии вызы- 
вает дополнительные издержки [2]. В связи с этим 
возникает необходимость расчета оптимальных 
значений диаметра газопровода Д для транспор- 
тировки заданного количества газа. 

Цель работы заключалась в том, чтобы на осно- 
ве предложенного уточненного метода расчета О ор 
определить погрешность известного приближен- 
ного алгоритма [3]. 

Предполагались известными следующие исход- 
ные данные: молекулярная масса перекачиваемого 
газа (метан, пентан и т.п.) М; его объемный расход 
при нормальных физических условиях Ѵ 0 , м 3 /с; 
длина горизонтального газопровода Ь, м; КПД на- 
соса г / „, передачи іу, двигателя г/з- Стоимости: 
электроэнергии 2е, за 1 кВт ч, амортизации 2а, и 
эксплуатации трубопровода 2г, руб. в год на 1 м 
длины и 1 м диаметра. 

В приближенной методике принималось, что ко- 
эффициент трения постоянен и равен Х р , и что поте- 
ри на местные сопротивления составляют к т (к т < 1) от 
потерь на трение [3]. Такое допущение, использован- 
ное в [3], объясняется, вероятно, тем, чтобы вычисле- 
ния не были итерационными. Как подчеркивается в 
[1], гидравлический расчет газопроводов следует вы- 
полнять с учетом зависимости коэффициентов гид- 
равлического сопротивления от режима движения 
газа (ламинарный, переходный, турбулентный). 


Расчет выполнялся для условной температуры 
Т. Принималось (с последующей проверкой), что 
падение давления в трубопроводе невелико, и 
среднее давление примерно равно атмосферному. 

Как известно [1], процесс оптимизации состоит 
из следующих операций: 1) составление целевой 
функции (ЦФ) и 2) отыскание значений парамет- 
ров, при которых ее значение будет экстремальным. 

Использовался экономический критерий опти- 
мальности, учитывающий приведенные затраты. 
Суммарная годовая стоимость газопровода как 
функция диаметра Д выражалась уравнением [3] 

Д = 5/ Д +(2а, + 2г х )Ш, 

Б=2 е „(1 + Д)(Я р /2Д-г г (4К/^) 2 , (1) 

где г р И- плотность газа и объемный расход газа; 2е„ 
- стоимость электроэнергии с учетом продолжи- 
тельности эксплуатации газопровода в течении года. 

С использованием аналитического условия ми- 
нимума ЦФ {дЕ ! /дГ)= 0, д 1 Е^д^ 2 > 0), в [3] получено 
выражение для оптимального значения диаметра 
газопровода 

О орр =\5В/(2а,+2 П )Ц' 6 . (2) 

Уточненный метод расчета основывался на том, 
что коэффициент сопротивления являлся функци- 
ей от П и вычислялся следующим образом [4]. 

Для труб круглого сечения при ламинарном ре- 
жиме течения газа (критерий Рейнольдса Ке<2300) - 

А„„=64/Ке, Ке = 4 рѴ ЦцП), (3) 

где /ц - коэффициент динамической вязкости. 

При турбулентном режиме течения (Ке > 2300) 
для шероховатых труб - 

\ь = ОД 1(Дкв / 2) + 68 / Ке) 0,25 , ( 4 ) 
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где к зт = 1,316^ - эквивалентная шероховатость, 
учитывающая как среднюю высоту выступов, так и 
их форму, расположение в плане; к„ - средняя вы- 
сота выступов шероховатости. 

Аналитический метод поиска экстремума уточ- 
ненной Цф, в которой коэффициент сопротивле- 
ния вычисляется по формулам (3) или (4), приво- 
дит к необходимости численного решения транс- 
цендентного уравнения. В данной работе для рас- 
чета минимума суммарной стоимости ( 1 ) исполь- 
зовался один из численных методов поиска [5]. 
При этом учитывалось свойство унимодальности 
ЦФ, т.е. наличие одного экстремума. Последова- 
тельно увеличивался диаметр трубы и вычислялось 
соответствующие значения ЦФ (1). Когда умень- 
шение 2" сменялось возрастанием, то цикл прекра- 
щался, и полученные значения П и Д принимались 
за искомые. Затем определялись относительные 
погрешности значений диаметров (е в ) и стоимости 
(в г ). Описанный алгоритм реализован в виде прог- 
раммного модуля в среде "ТигЬо РазсаГ. 

Расчет приведенного в [3] примера показал сле- 
дующее. При транспортировке метана, имеющего 
температуру 7=303 К, объемный расход К 0 =2,16 м 3 /с, 
на расстояние 4 км приближенное значение опти- 
мального диаметра Д )д , по формуле (2) равно 0,663 м. 
Расчет по уточненной методике определил значение 

= 0,595 м. Погрешность при проведении сравни- 
тельного анализа составила 11,4 %. Среди труб, вы- 
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пускаемых промышленностью, ближай ш ими к по- 
лученным значениям будут трубы с наружными диа- 
метрами 0,630 и 0,720м [1]. Реальное отличие состав- 
ляет уже 14,3 %. 

Суммарные затраты по уточненной методике на 
2,5 % ниже. С увеличением объемного расхода газа 
погрешность приближенного метода возрастает 
(таблица). 

Таблица. Результаты расчета оптимального значения диа- 
метра газопровода по приближенной и точной 
методиках 


К. м 7с 

1,85 

2,0 

2,15 

2,30 

2,45 

°о № . М 

0,614 

0,639 

0,663 

0,686 

0,709 

С„„, м 

0,554 

0,575 

0,595 

0,616 

0,635 

Е в , % 

10,8 

11,1 

11,4 

11,5 

11,7 

е„ % 

2,5 

2,6 

2,7 

2,7 

2,8 


Таким образом, установлено, что приближен- 
ная методика дает завышенные значения опти- 
мального диаметра и суммарных затрат газопрово- 
да. Поэтому корректнее (даже в учебных целях) 
проводить расчет по уточненной методике. 

В курсе "Надежность и оптимизация систем 
теплоэнергоснабжения промышленных предприя- 
тий", читаемом на кафедре "Теоретической и про- 
мышленной теплотехники" использован предло- 
женный алгоритм. 
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